Discovery and characterization of novel DNA viruses in Apis mellifera: expanding the honey bee virome through metagenomic analysis.

To date, many viruses have been discovered to infect honey bees. In this study, we used high-throughput sequencing to expand the known virome of the honey bee, Apis mellifera, by identifying several novel DNA viruses. While the majority of previously identified bee viruses are RNA, our study reveals nine new genomes from the Parvoviridae family, tentatively named Bee densoviruses 1 to 9. In addition, we characterized a large DNA virus, Apis mellifera filamentous-like virus (AmFLV), which shares limited protein identities with the known Apis mellifera filamentous virus. The complete sequence of AmFLV, obtained by a combination of laboratory techniques and bioinformatics, spans 152,678 bp. Linear dsDNA genome encodes for 112 proteins, of which 49 are annotated. Another large virus we discovered is Apis mellifera nudivirus, which belongs to a group of Alphanudivirus. The virus has a length of 129,467 bp and a circular dsDNA genome, and has 106 protein encoding genes. The virus contains most of the core genes of the family Nudiviridae. This research demonstrates the effectiveness of viral binning in identifying viruses in honey bee virology, showcasing its initial application in this field.IMPORTANCEHoney bees contribute significantly to food security by providing pollination services. Understanding the virome of honey bees is crucial for the health and conservation of bee populations and also for the stability of the ecosystems and economies for which they are indispensable. This study unveils previously unknown DNA viruses in the honey bee virome, expanding our knowledge of potential threats to bee health. The use of the viral binning approach we employed in this study offers a promising method to uncovering and understanding the vast viral diversity in these essential pollinators.

Viral communities in millipede guts: Insights into the diversity and potential role in modulating the microbiome.

Millipedes are important detritivores harbouring a diverse microbiome. Previous research focused on bacterial and archaeal diversity, while the virome remained neglected. We elucidated the DNA and RNA viral diversity in the hindguts of two model millipede species with distinct microbiomes: the tropical Epibolus pulchripes (methanogenic, dominated by Bacillota) and the temperate Glomeris connexa (non-methanogenic, dominated by Pseudomonadota). Based on metagenomic and metatranscriptomic assembled viral genomes, the viral communities differed markedly and preferentially infected the most abundant prokaryotic taxa. The majority of DNA viruses were Caudoviricetes (dsDNA), Cirlivirales (ssDNA) and Microviridae (ssDNA), while RNA viruses consisted of Leviviricetes (ssRNA), Potyviridae (ssRNA) and Eukaryotic viruses. A high abundance of subtypes I-C, I-B and II-C CRISPR-Cas systems was found, primarily from Pseudomonadota, Bacteroidota and Bacillota. In addition, auxiliary metabolic genes that modulate chitin degradation, vitamins and amino acid biosynthesis and sulphur metabolism were also detected. Lastly, we found low virus-to-microbe-ratios and a prevalence of lysogenic viruses, supporting a Piggyback-the-Winner dynamic in both hosts.

Unveiling CRESS DNA Virus Diversity in Oysters by Virome.

Oysters that filter feed can accumulate numerous pathogens, including viruses, which can serve as a valuable viral repository. As oyster farming becomes more prevalent, concerns are mounting about diseases that can harm both cultivated and wild oysters. Unfortunately, there is a lack of research on the viruses and other factors that can cause illness in shellfish. This means that it is harder to find ways to prevent these diseases and protect the oysters. This is part of a previously started project, the Dataset of Oyster Virome, in which we further study 30 almost complete genomes of oyster-associated CRESS DNA viruses. The replication-associated proteins and capsid proteins found in CRESS DNA viruses display varying evolutionary rates and frequently undergo recombination. Additionally, some CRESS DNA viruses have the capability for cross-species transmission. A plethora of unclassified CRESS DNA viruses are detectable in transcriptome libraries, exhibiting higher levels of transcriptional activity than those found in metagenome libraries. The study significantly enhances our understanding of the diversity of oyster-associated CRESS DNA viruses, emphasizing the widespread presence of CRESS DNA viruses in the natural environment and the substantial portion of CRESS DNA viruses that remain unidentified. This study's findings provide a basis for further research on the biological and ecological roles of viruses in oysters and their environment.

Exploring the gut DNA virome in fecal immunochemical test stool samples reveals associations with lifestyle in a large population-based study.

Stool samples for fecal immunochemical tests (FIT) are collected in large numbers worldwide as part of colorectal cancer screening programs. Employing FIT samples from 1034 CRCbiome participants, recruited from a Norwegian colorectal cancer screening study, we identify, annotate and characterize more than 18000 DNA viruses, using shotgun metagenome sequencing. Only six percent of them are assigned to a known taxonomic family, with Microviridae being the most prevalent viral family. Linking individual profiles to comprehensive lifestyle and demographic data shows 17/25 of the variables to be associated with the gut virome. Physical activity, smoking, and dietary fiber consumption exhibit strong and consistent associations with both diversity and relative abundance of individual viruses, as well as with enrichment for auxiliary metabolic genes. We demonstrate the suitability of FIT samples for virome analysis, opening an opportunity for large-scale studies of this enigmatic part of the gut microbiome. The diverse viral populations and their connections to the individual lifestyle uncovered herein paves the way for further exploration of the role of the gut virome in health and disease.

Distribution patterns and functional diversity of DNA viruses determined by ecological niches in huge river ecosystems.

While the vast number of DNA and RNA viruses participate in biogeochemical cycles in natural systems, little is known about virome in river ecosystems. Here, we analyzed the DNA viral composition and its metabolic potential in the Yangtze River, including freshwater (FW) and freshwater sediments (FWS). A total of 1237 river-derived virus contigs (RVCs) were obtained following de novo assembly from 62 metagenomics. We found that the viral diversity is significantly positively correlated longitudinally. Moreover, FW exhibited a greater viral variety and significantly different composition than FWS. The viral co-occurrence network suggested that positive correlations predominate between RVCs. Lastly, 1657 viral functions were predicted by gene ontology. Notably, 96 of 150 RVCs with higher weights identified by random-forest classier were more abundant in FW, which most engage organic cyclic compound metabolic processes and hydrolase activity. Together, this study highlights the previously unrecognized viruses and the importance of their distributions and functions in major river systems.

In vitro characterization of cell-fusing agent virus DNA forms in Aedes aegypti mosquitoes.

How non-retroviral endogenous viral elements (EVEs) are established is a long-standing question. Viral DNA (vDNA) forms of RNA viruses are likely to be EVE precursors. Cell-fusing agent virus (CFAV) is a major insect-specific virus (ISV) in the Aedes aegypti mosquitoes and one of the few existing non-retroviral RNA viruses found as EVEs. We characterized CFAV-derived vDNA in the cell line to understand the mechanism of why current viruses are rarely endogenized. vDNA production was affected by cell culture media independent of CFAV replication. vDNAs that correspond to different regions covering the entire viral genome were detected, implying multiple initiation sites exist. A considerable proportion of vDNAs corresponded to ssDNA. Higher vDNA copies were detected in the cytoplasm than the nucleus. Our findings provide valuable insights into the intracellular characteristics of ISV-derived vDNAs, which will aid in understanding the underlying mechanisms of non-retroviral EVE formation.

Environmental viromes reveal the global distribution signatures of deep-sea DNA viruses.

INTRODUCTION: Viruses are abundant and ecologically significant in marine ecosystems. However, the virome of deep-sea sediments is not extensively investigated. OBJECTIVES: To explore the distribution pattern of deep-sea viruses on a global scale, the viromes of DNA viruses isolated from 138 sediments of 5 deep-sea ecosystems were characterized. METHODS: The viral particles were purified from each sediment sample. Then the viral DNAs were extracted and subjected to viral metagenomic analysis. RESULTS: Here, we constructed a global deep-sea environmental virome dataset by analyzing the viral DNA of 138 sediment samples. A total of 347,737 viral operational taxonomic units (vOTUs) were identified, of which 84.94% were hitherto unknown, indicating that deep sea was a reservoir of novel DNA viruses. Furthermore, circular viral genome analysis revealed 98,581 complete genomes. The classified vOTUs included eukaryotic (44.55%) and prokaryotic (25.75%) viruses, and were taxonomically assigned to 63 viral families. The composition and abundance of the deep-sea sediment viromes were dependent on the deep-sea ecosystem as opposed to geographical region. Further analysis revealed that the viral community differentiation in different deep-sea ecosystems was driven by the virus-mediated energy metabolism. CONCLUSION: Our findings showed that deep-sea ecosystems are a reservoir of novel DNA viruses and the viral community is shaped by the environmental characteristics of deep-sea ecosystems, thus presenting critical information for determining the ecological significance of viruses in global deep-sea ecosystems.

Integrase-associated niche differentiation of endogenous large DNA viruses in crustaceans.

IMPORTANCE: Crustacean genomes harbor sequences originating from a family of large DNA viruses called nimaviruses, but it is unclear why they are present. We show that endogenous nimaviruses selectively insert into repetitive sequences within the host genome, and this insertion specificity was correlated with different types of integrases, which are DNA recombination enzymes encoded by the nimaviruses themselves. This suggests that endogenous nimaviruses have colonized various genomic niches through the acquisition of integrases with different insertion specificities. Our results point to a novel survival strategy of endogenous large DNA viruses colonizing the host genomes. These findings may clarify the evolution and spread of nimaviruses in crustaceans and lead to measures to control and prevent the spread of pathogenic nimaviruses in aquaculture settings.

Circular replication-associated protein-encoding single-stranded DNA virus with risk of spillover is widely prevalent in domestic animals in China.

Circular replication-associated protein (Rep)-encoding single-stranded (CRESS) DNA viruses are highly diverse and have a broad range of hosts. In this study, we report the detection of Bo-Circo-like virus AH20-1 in the feces of diarrheal cattle. The virus has a circular genome of 3,912 nucleotides, three major putative open reading frames, and encodes a Rep gene of 310 amino acids. We found that the virus is closely related to the Bo-Circo-like virus CH strain, which belongs to the novel Kirkoviridae family. Furthermore, we conducted a nationwide surveillance program and found that the virus is prevalent in China (23.6%, 205/868), with the BCLa subtype being the predominant strain. Our findings suggest that the virus can infect sheep, highlighting the potential for cross-species transmission. Our pressure analysis indicates that the CRESS-DNA Kirkoviridae family has broad host adaptation, and that selection pressure played an important role in the evolution of its Rep genes. Our study underscores the need for continued epidemiological surveillance of this virus due to its widespread prevalence in our ruminant population and potential for cross-species transmission.

Ostreid herpesvirus 1 latent infection and reactivation in adult Pacific oysters, Crassostrea gigas.

Ostreid herpesvirus 1 (OsHV-1) is one of the most economically important pathogens of Pacific oysters. Understanding the pathogenesis of this virus is critical to developing tools to control outbreaks on shellfish farms. OsHV-1 is genetically related to vertebrate herpesviruses, which have a lytic and a latent stage, with the latent stage capable of being reactivated to the lytic stage. Here, OsHV-1 latency in Pacific oysters was investigated in experimentally and naturally infected oysters. Lytic infection in one-year-old oysters injected with the Tomales Bay strain of OsHV-1 was detectable between 1 and 4 days post-injection (dpi) but was not detectable after 5 dpi. The injected oysters shed 1 × 102 to 1 × 104 DNA copies/ml into the water during the 4-day acute phase. Lytic shedding was not detectable in two-year-old oysters injected similarly with the same strain of OsHV-1; however, the OsHV-1 genome was detectable by qPCR in the adductor muscle, gill, mantle, and hemocytes within the first 3 dpi, after which it became undetectable. No OsHV-1 was detectable in the adductor muscle, gill, or mantle from experimentally infected oysters on days 15 and 21 post-injection or from oysters sampled 9 months after surviving an OsHV-1 mortality event; however, OsHV-1 DNA could be detected in hemocytes of both experimentally infected oysters at 21 dpi and naturally infected oysters using nested PCR. In addition, lytic viral gene transcription was detectable in hemocytes of experimentally infected oysters between 1 and 21 dpi and in hemocytes of naturally infected oysters. Furthermore, OsHV-1 reactivation from latency was induced in experimentally infected oysters at 21 dpi and in naturally infected oysters 12 months after an OsHV-1 outbreak.

Identification of small circular DNA viruses in coyote fecal samples from Arizona (USA).

Coyotes (Canis latrans) have a broad geographic distribution across North and Central America. Despite their widespread presence in urban environments in the USA, there is limited information regarding viruses associated with coyotes in the USA and in particular the state of Arizona. To explore viruses associated with coyotes, particularly small DNA viruses, 44 scat samples were collected (April-June 2021 and November 2021-January 2022) along the Salt River near Phoenix, Arizona (USA), along 43 transects (500 m). From these samples, we identified 11 viral genomes: two novel circoviruses, six unclassified cressdnaviruses, and two anelloviruses. One of the circoviruses is most closely related to a circovirus sequence identified from an aerosolized dust sample in Arizona, USA. The second circovirus is most closely related to a rodent-associated circovirus and canine circovirus. Of the unclassified cressdnaviruses, three encode replication-associated proteins that are similar to those found in protists (Histomonas meleagridis and Monocercomonoides exilis), implying an evolutionary relationship with or a connection to similar unidentified protist hosts. The two anelloviruses are most closely related to those found in rodents, and this suggests a diet-related identification.

Novel plant-associated genomoviruses from the Brazilian Cerrado biome.

The discovery rate of new plant viruses has increased due to studies involving high-throughput sequencing (HTS), particularly for single-stranded DNA viruses of the family Genomoviridae. We carried out an HTS-based survey of genomoviruses in a wide range of native and exotic trees grown in the Brazilian Cerrado biome, and the complete genome sequences of two novel members of the family Genomoviridae from two distinct genera were determined. Specific primers were designed to detect these genomoviruses in individual samples. A new gemykolovirus (Tecoma stans associated gemykolovirus) was detected in Tecoma stans, and a new gemykibivirus (Ouratea duparquetiana associated gemykibivirus) was detected in Ouratea duparquetiana. A gemykrogvirus related to Gila monster associated gemykrogvirus (80% pairwise identity) was also detected in foliar samples of Trembleya parviflora. Our pilot study paves the way for a better characterization of this diverse collection of genomoviruses as well as their interactions with the associated tree species.

Genomic adaptation of giant viruses in polar oceans.

Despite being perennially frigid, polar oceans form an ecosystem hosting high and unique biodiversity. Various organisms show different adaptive strategies in this habitat, but how viruses adapt to this environment is largely unknown. Viruses of phyla Nucleocytoviricota and Mirusviricota are groups of eukaryote-infecting large and giant DNA viruses with genomes encoding a variety of functions. Here, by leveraging the Global Ocean Eukaryotic Viral database, we investigate the biogeography and functional repertoire of these viruses at a global scale. We first confirm the existence of an ecological barrier that clearly separates polar and nonpolar viral communities, and then demonstrate that temperature drives dramatic changes in the virus-host network at the polar-nonpolar boundary. Ancestral niche reconstruction suggests that adaptation of these viruses to polar conditions has occurred repeatedly over the course of evolution, with polar-adapted viruses in the modern ocean being scattered across their phylogeny. Numerous viral genes are specifically associated with polar adaptation, although most of their homologues are not identified as polar-adaptive genes in eukaryotes. These results suggest that giant viruses adapt to cold environments by changing their functional repertoire, and this viral evolutionary strategy is distinct from the polar adaptation strategy of their hosts.

Benchmarking of virome metagenomic analysis approaches using a large, 60+ members, viral synthetic community.

IMPORTANCE: We report here efforts to benchmark performance of two widespread approaches for virome analysis, which target either virion-associated nucleic acids (VANA) or highly purified double-stranded RNAs (dsRNAs). This was achieved using synthetic communities of varying complexity levels, up to a highly complex community of 72 viral agents (115 viral molecules) comprising isolates from 21 families and 61 genera of plant viruses. The results obtained confirm that the dsRNA-based approach provides a more complete representation of the RNA virome, in particular, for high complexity ones. However, for viromes of low to medium complexity, VANA appears a reasonable alternative and would be the preferred choice if analysis of DNA viruses is of importance. Several parameters impacting performance were identified as well as a direct relationship between the completeness of virome description and sample sequencing depth. The strategy, results, and tools used here should prove useful in a range of virome analysis efforts.

Taxonomic update for giant viruses in the order Imitervirales (phylum Nucleocytoviricota).

Large DNA viruses in the phylum Nucleocytoviricota, sometimes referred to as "giant viruses" owing to their large genomes and virions, have been the subject of burgeoning interest over the last decade. Here, we describe recently adopted taxonomic updates for giant viruses within the order Imitervirales. The families Allomimiviridae, Mesomimiviridae, and Schizomimiviridae have been created to accommodate the increasing diversity of mimivirus relatives that have sometimes been referred to in the literature as "extended Mimiviridae". In addition, the subfamilies Aliimimivirinae, Megamimivirinae, and Klosneuvirinae have been established to refer to subgroups of the Mimiviridae. Binomial names have also been adopted for all recognized species in the order. For example, Acanthamoeba polyphaga mimivirus is now classified in the species Mimivirus bradfordmassiliense.

Global diversity and biogeography of DNA viral communities in activated sludge systems.

BACKGROUND: Activated sludge (AS) systems in wastewater treatment plants (WWTPs) harbor enormous viruses that regulate microbial metabolism and nutrient cycling, significantly influencing the stability of AS systems. However, our knowledge about the diversity of viral taxonomic groups and functional traits in global AS systems is still limited. To address this gap, we investigated the global diversity and biogeography of DNA viral communities in AS systems using 85,114 viral operational taxonomic units (vOTUs) recovered from 144 AS samples collected across 54 WWTPs from 13 different countries. RESULTS: AS viral communities and their functional traits exhibited distance-decay relationship (DDR) at the global scale and latitudinal diversity gradient (LDG) from equator to mid-latitude. Furthermore, it was observed that AS viral community and functional gene structures were largely driven by the geographic factors and wastewater types, of which the geographic factors were more important. Carrying and disseminating auxiliary metabolic genes (AMGs) associated with the degradation of polysaccharides, sulfate reduction, denitrification, and organic phosphoester hydrolysis, as well as the lysis of crucial functional microbes that govern biogeochemical cycles were two major ways by which viruses could regulate AS functions. It was worth noting that our study revealed a high abundance of antibiotic resistance genes (ARGs) in viral genomes, suggesting that viruses were key reservoirs of ARGs in AS systems. CONCLUSIONS: Our results demonstrated the highly diverse taxonomic groups and functional traits of viruses in AS systems. Viral lysis of host microbes and virus-mediated HGT can regulate the biogeochemical and nutrient cycles, thus affecting the performance of AS systems. These findings provide important insights into the viral diversity, function, and ecology in AS systems on a global scale. Video Abstract.

Diversity and function of mountain and polar supraglacial DNA viruses.

Mountain and polar glaciers cover 10% of the Earth's surface and are typically extreme environments that challenge life of all forms. Viruses are abundant and active in supraglacial ecosystems and play a crucial role in controlling the supraglacial microbial communities. However, our understanding of virus ecology on glacier surfaces and their potential impacts on downstream ecosystems remains limited. Here, we present the supraglacial virus genome (SgVG) catalog, a 15-fold expanded genomic inventory of 10,840 DNA-virus species from 38 mountain and polar glaciers, spanning habitats such as snow, ice, meltwater, and cryoconite. Supraglacial DNA-viruses were highly specific compared to viruses in other ecosystems yet exhibited low public health risks. Supraglacial viral communities were primarily constrained by habitat, with cryoconite displaying the highest viral activity levels. We observed a prevalence of lytic viruses in all habitats, especially in cryoconite, but a high level of lysogenic viruses in snow and ice. Additionally, we found that supraglacial viruses could be linked to ∼83% of obtained prokaryotic phyla/classes and possessed the genetic potential to promote metabolism and increase cold adaptation, cell mobility, and phenolic carbon use of hosts in hostile environmental conditions using diverse auxiliary metabolic genes. Our results provide the first systematic characterization of the diversity, function, and public health risks evaluation of mountain and polar supraglacial DNA viruses. This understanding of glacial viruses is crucial for function assessments and ecological modeling of glacier ecosystems, especially for the Tibetan Plateau's Mountain glaciers, which support ∼20% of the human populations on Earth.

Pangenomic Analysis of Nucleo-Cytoplasmic Large DNA Viruses. I: The Phylogenetic Distribution of Conserved Oxygen-Dependent Enzymes Reveals a Capture-Gene Process.

The Nucleo-Cytoplasmic Large DNA Viruses (NCLDVs) infect a wide range of eukaryotic species, including amoeba, algae, fish, amphibia, arthropods, birds, and mammals. This group of viruses has linear or circular double-stranded DNA genomes whose size spans approximately one order of magnitude, from 100 to 2500 kbp. The ultimate origin of this peculiar group of viruses remains an open issue. Some have argued that NCLDVs' origin may lie in a bacteriophage ancestor that increased its genome size by subsequent recruitment of eukaryotic and bacterial genes. Others have suggested that NCLDVs families originated from cells that underwent an irreversible process of genome reduction. However, the hypothesis that a number of NCLDVs sequences have been recruited from the host genomes has been largely ignored. In the present work, we have performed pangenomic analyses of each of the seven known NCLDVs families. We show that these families' core- and shell genes have cellular homologs, supporting possible escaping-gene events as part of its evolution. Furthermore, the detection of sequences that belong to two protein families (small chain ribonucleotide reductase and Erv1/Air) and to one superfamily [2OG-Fe(II) oxygenases] that are for distribution in all NCLDVs core and shell clusters encoding for oxygen-dependent enzymes suggests that the highly conserved core these viruses originated after the Proterozoic Great Oxidation Event that transformed the terrestrial atmosphere 2.4-2.3 Ga ago.

Characterization, Directed Evolution, and Targeting of DNA Virus-Encoded RNA Capping Enzymes Using Phenotypic Yeast Platforms.

The constant and the sudden emergence of zoonotic human and animal viruses is a significant threat to human health, the world economy, and the world food supply. This has necessitated the development of broad-spectrum therapeutic strategies to combat these emerging pathogens. Mechanisms that are essential for viral replication and propagation have been successfully targeted in the past to develop broad-spectrum therapeutics that can be readily repurposed to combat new zoonotic pathogens. Because of the importance of viral RNA capping enzymes to viral replication and pathogenesis, as well as their presence in both DNA and RNA viruses, these viral proteins have been a long-standing therapeutic target. Here, we use genome sequencing information and yeast-based platforms (YeRC0M) to identify, characterize, and target viral genome-encoded essential RNA capping enzymes from emerging strains of DNA viruses, i.e., Monkeypox virus and African Swine Fever Virus, which are a significant threat to human and domestic animal health. We first identified and biochemically characterized these viral RNA capping enzymes and their necessary protein domains. We observed significant differences in functional protein domains and organization for RNA capping enzymes from emerging DNA viruses in comparison to emerging RNA viruses. We also observed several differences in the biochemical properties of these viral RNA capping enzymes using our phenotypic yeast-based approaches (YeRC0M) as compared to the previous in vitro studies. Further, using directed evolution, we were able to identify inactivation and attenuation mutations in these essential viral RNA capping enzymes; these data could have implications on virus biocontainment as well as live attenuated vaccine development. We also developed methods that would facilitate high-throughput phenotypic screening to identify broad-spectrum inhibitors that selectively target viral RNA capping enzymes over host RNA capping enzymes. As demonstrated here, our approaches to identify, characterize, and target viral genome-encoded essential RNA capping enzymes are highly modular and can be readily adapted for targeting emerging viral pathogens as well as their variants that emerge in the future.

Rapid and sensitive single-sample viral metagenomics using Nanopore Flongle sequencing.

The ability of viral metagenomic Next-Generation Sequencing (mNGS) to unbiasedly detect nucleic acids in a clinical sample is a powerful tool for advanced diagnosis of viral infections. When clinical symptoms do not provide a clear differential diagnosis, extensive laboratory testing with virus-specific PCR and serology can be replaced by a single viral mNGS analysis. However, widespread diagnostic use of viral mNGS is thus far limited by long sample-to-result times, as most protocols rely on Illumina sequencing, which provides high and accurate sequencing output but is time-consuming and expensive. Here, we describe the development of an mNGS protocol based on the more cost-effective Nanopore Flongle sequencing with decreased turnaround time and lower, yet sufficient sequencing output to provide sensitive virus detection. Sample preparation (6 h) and sequencing (2 h) times are substantially reduced compared to Illumina mNGS and allow detection of DNA/RNA viruses at low input (up to 33-38 cycle threshold of specific qPCR). Although Flongles yield lower sequencing output, direct comparison with Illumina mNGS on diverse clinical samples showed similar results. Collectively, the novel Nanopore mNGS approach is specifically tailored for use in clinical diagnostics and provides a rapid and cost-effective mNGS strategy for individual testing of severe cases.